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ABSTRACT 

We present observations with the IRAM Plateau de Bure Interferometer of three QSOs at z>5 aimed at detecting molecular gas in 
their host galaxies as traced by CO transitions. CO (5-4) is detected in SDSS J033829.31+002156.3 at z=5.0267, placing it amongst 
the most distant sources detected in CO. The CO emission is unresolved with a beam size of ~ 1", implying that the molecular gas 
is contained within a compact region, less than ~ 3 kpc in radius. We infer an upper limit on the dynamical mass of the CO emitting 
region of ~ 3 x 10 10 Mq/ sin(i) 2 . The comparison with the Black Hole mass inferred from near-IR data suggests that the BH-to- 
bulge mass ratio in this galaxy is significantly higher than in local galaxies. From the CO luminosity we infer a mass reservoir of 
molecular gas as high as M(H 2 ) = 2.4 x 10'°Mq, implying that the molecular gas accounts for a significant fraction of the dynamical 
mass. When compared to the star formation rate derived from the far-IR luminosity, we infer a very short gas exhaustion timescale 
(~ 10 7 years), comparable to the dynamical timescale. CO is not detected in the other two QSOs (SDSS J083643. 85+005453.3 and 
SDSS J163033. 90+401209. 6) and upper limits are given for their molecular gas content. When combined with CO observations of 
other type 1 AGNs, spanning a wide redshift range (0<z<6.4), we find that the host galaxy CO luminosity (hence molecular gas con- 
tent) and the AGN optical luminosity (hence BH accretion rate) are correlated, but the relation is not linear: L£, oc [/U_a(4400A)] 



Moreover, at high redshifts (and especially at z>5) the CO luminosity appears to saturate. We discuss the implications of these findings 
in terms of black hole-galaxy co-evolution. 

Key words. Galaxies: high redshift - Galaxies: ISM - quasars: general - Infrared: galaxies - Submillimeter 

1. Introduction ISolomon et al.L fl997l) : hence strong far-IR emission may be a 

_ , . . , . , . . .... , pre-requisite for CO detection. However, at z>5 the ISM is ex- 

The detection of carbon monoxide (CO) emission in high red- d tQ und s met allicity and dust evolution, which 

shift galaxies provides a crucial tool for investigating the early cause M h . z , axies tQ deyiate from ±e ^ CQ _ mR fe _ 

epochs of galaxy formation (see |Solomon & Vanden Bern* ^on. Another caveat is that present mm/submm detections are 

|2005l for a review). Indeed, CO emission is a proxy for the doge tQ ^ sensitivit limit of current cameras; hence even a 

molecular gas content, the reservoir for star formation. The smaU scatter fa ±e CQ/FlR mio lead tQ a non . detection in 
CO line profile also provides information on galaxy kinematics, 

from which constraints on the dynamical mass can be inferred. FIR 

Currently only 9 galaxies have been detected in CO at z>4, An additional issue affecting the detection of CO in high- 
only three of which are at z>5, and all of them host po werful z Q SOs 18 the accuracy of the redshift. Indeed, at z>4 the 
QSOs or radio galaxies ISotorrTon & Vanden Boul EPOS). With emission lines typically observed in the optical band are ei- 
the ex ception of the radio-galaxy TN J0924-2201 klamer et all mer strongly blueshifted with respect to the systemic velocity 
2005), all of these high-z CO detections were obtained in galax- of the host gal axy (such as CIV at 1549 A and SilV at 1400A, 



ies pre-selected amongst luminous far-IR sources (as inferred |Richards et aU |2002|) or, in the case of Lya, strongly affected 

from mm/submm continuum observations). However, such a by intergalactic gas absorption (Lya Forest). In these cases the 

selection criterion may prevent us from identifying evolution- redshift deviations from the systemic velocity can be as large 

ary effects in terms of molecular gas and dust content in as several thousand km s 1 . Until recently, millimetre receivers 

high-z galaxies. Indeed, in local gal axies, CO and FIR lumi- had bandwidths of 0.5 GHz (covering 1500 km s 1 at the best), 

nosities are known to correlate (e.g. lYoung & Scovilld Il99lt limiting the efficiency of the search for CO in sources with such 

uncertain redshift estimates. This has changed very recently with 

Send offprint requests to: R. Maiolino the implementation of new receivers having a much wider band- 
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width (e.g. the new receivers at the Plateau de Bure have a band- 
width of 4 GHz). Additionally, one can observe lower ionization 
lines, such as Mgll at 2798 A and CIII] at 1909 A, which pro- 
vide a better redshift estimate, since they are generally shifted 
by only a few hundred km s _1 with respect to the systemic ve- 
locity of the host galaxy. At z>4 such low-ionization lines are 
shifted into the near-IR, and near-IR spectroscopic campaigns 
have recently provided accurate redshifts for a number of high-z 
QSOs. 

With the goal of increasing the number of CO detections 
at z>5, and removing any bias towards FIR-luminous sources, 
we observed three high-z QSOs with the Plateau de Bure 
Interferometer (PdBI). Here we report the detection of CO emis- 
sion in one source, and upper limits for the other two. We adopt 
the following cosmological parameters: Hp — 7 1 km s" 1 Mpc -1 , 
Q A = 0.73 and Q m = 0.27 dSpergel et all 120031) . 



2. Sample selection and observations 

The three QSOs were selected from the SDSS catalog to be at 
z>5 and observable from the IRAM PdBI. For all, the redshift 
had been re-determined with Mgll or CIII] near-IR spectroscopy, 
as listed in Tabl^T] These QSOs were selected regardless of their 
FIR luminosity, as inferred from previous submm/mm bolomet- 
ric observations: two of the sources have not been detected in 
continuum at 850/rni nor at 1.2mm (Tabl£D», while J0338+0021 
has a detection both at 85 0/jm and at 1 .2mm, fr om which Lfir = 

I. 5 x 10 13 L© is inferred (Priddev et al., 2003). Finally, we note 
that J0836+0054 is a radio loud QSO dPetric et al.Ll2003l) . 

Observations in configuration D were performed between 
November 2005 and July 2006 with the IRAM PdBI six ele- 
ments array. The old generation 3mm receivers were tuned in 
single sideband to the frequency of the redshifted CO(5-4) or 
CO(6-5) line, depending on the specific redshift of each source 
(Table [T). The beam size in D configuration at such frequen- 
cies is typically ~ 5". The on-source integration times were 

II. 2 hours for J0338+0021, 10.5 hours for J0836+0054 and 12.6 
hours for J1630+4012. 

Following the CO detection in J0338+0021 (see next sec- 
tion), this source was also observed in configuration A (with the 
new PdBI receivers). Observations were obtained in January and 
February 2007, for a total of 5.9 hours on source. The resulting 
synthesized beam size is 1.59" x 0.76" (position angle: 28°). 

The resulting lcr sensitivities are 0.23, 0.34 and 0.26 
mJy/beam in channels of 400 km s _1 for J0338+0021, 
J0836+0054 and J1630+4012, respectively. 



3. Results 

CO(5^1) is detected in J0338+0021 with a significance of 
8cr, at a frequency of 95.6191 GHz. This is amongst the 
most distant CO detections obtained so far, together with 
SDSS Jl 14816.64+525150.3 at z=6.4, TN J 0924-2201 at z=5.2 
and SDSS J09272 1.82+200 123 .7 at z=5.7 7 JWaite7et^l200l 
iBertoldi et all 120031: iKlamer et all 120051: ICarilli et all l2007h . 
Fig- [U shows the spectrum and Figf5]presents the integrated in- 
tensity map. The rms per channel of the spectrum in Fig. Q] 
is 0.55 mJy/beam. We note that the individual spectra taken 
in A and D configuration are consistent with each other. The 
redshift inferred from the CO line is z=5.0267±0.0003, i.e. 
consistent with the Mgll redshift (Tabl^B within the uncer- 
tainty of ~20 km 
sition of the CO 



s 1 on the line center. The absolute po- 
source is at RA(J200)=03:38:29.32 and 



DEC(J2000)=00:2 1:56.1 (accuracy < 0.1"), which is consistent 
with the optical position. The source is spatially unresolved, im- 
plying a radius smaller than ~2.5 kpc (along the beam minor 
axis). 

In TableEpwe also report the CO luminosity L^, defined as in 
ISolomon et alj dl997l) . The millimetric continuum at the location 
ofJ0338+0021 is undetected (S v [95.6GHz] = -18±85/Jy); the 
upper limit is consistent with the extrapolation of the detections 
at higher freq uencies (using the g rey body fitting curves with 
T-50-65K in PriddexeLalillQOl. 



Both J0836+0054 and J1630+4012 (i.e. the two FIR-faint 
QSOs) were undetected in CO, and Table Q] gives the inferred 
upper limits on the CO intensity and luminosity (assuming a line 
width of 400 km s" 1 ). 



Redshift 
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Fig. 1. CO(5^1) spectrum of SDSSJ0338+0021 (sum of A and 
D conf.) rebinned to 10 MHz (31.4 km s~'). The red line shows 
a single gaussian fit to the line. Velocities are relative to the Mgll 
redshift. 
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Fig. 2. CO(5^) cleaned map of SDSSJ0338+0021 (from both 
configurations A and D) obtained by integrating in velocity from 
-300 to +230 km s _1 with contours in steps of and starting at 
lcr = 0.14 Jy beam -1 km s _1 . 
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4. Discussion 

4.1. Gas and dynamical mass 

The CO luminosity can be used to infer the molecular 
gas content. Most authors adopt the conversion factor a = 
0.8 M Q (K km s -1 pc 2 )~' between CO(l-O) line luminosity 
and M(H2), as inferred for nearby s tarbursts with moderate CO 
excitation and in virial equilibri um dDownes & Solomon! 1 1998t 
ISolomon & Vanden Bou t, 2005). Such a conversion factor may 
not be appropriate for high-z QSOs, which often show indica- 
tions of high gas excitation dBertoldi et all l2003t IWeiB et all 
I2007I) . We do not have information on the CO transitions lower 
than (5-4) in J0338+0021, and hence we cannot constrain the 
gas excitation. Given these uncertainties, we assume the same 
conversion factor as for local starbursts, and we also assume con- 
stant brightness temperature (the optically thick case) from J=l 
to J=5, i.e. Lp (l -0) = Lq (5 -4). Under these assumptions we 
infer a molecular gas mass of ~ 2.2 x 10 10 Mq. Note that this is 
probably a lower limit on the H2 mass, since the transition (5—4) 
is likely subthermal and the conversion factor probably higher, 
as inferred in other starburst galaxies and powerful sources (e.g. 
iBavet et alll20MlBradfordet^l200llWeiB et aUl2007t) . 

The line width and the upper limit on the CO size of 
J0338+0021 allow u s to infer an upper l imit on the dynam- 
ical mass. Following iTacconi et al.l (|2006) we derive v c sin(i) 
(where v c is the circular velocity at the outer CO radius, and 
"i" the inclination angle of the gaseous disk) by dividing the 
CO line FWHM by 2.4. We obtain M dyn = Rv 2 /G < 3.2 x 
10 10 M©/sin 2 (i), where we have assumed an upper limit for the 
size of the CO source of 1" (~ average of our beam sizes), hence 
R<3.2 kpc. The main uncertainty of the dynamical mass upper 
li mit is due to the unkno wn inclination angle "i". As discussed 
in ICarilli & Wangl ((2006), type 1 AGNs may be biased against 
edge-on host disks, because in such cases the nucleus should 
be obscured. Such a bias is inferred from the finding that the 
CO emission in (type 1) QSOs is systematically narrower than 
in SMGs at similar redshifts. However, we note that in the spe- 
cific case of J0338+0021 the CO FWHM is amongst the largest 
ever observed in QSO s (whose median FWHM is 300 km s -1 , 
ICarilli & Wangl 120061) . and similar to the median FWHM ob- 
served in SMGs (500 km s" 1 ), suggesting that J0338+0021 is 
likely observed at high inclination. 

Based on these results we conclude that the molecular gas 
mass accounts for a substantial fraction of the dynamical mass. 
More specifically, if the system is nearly edge-on (i ~ 90°) then 
the molecular gas mass accounts for more than 70% of the dy- 
namical mass. Even if the system has an inclination of i = 30° 
the molecular gas mass still accounts for more than 20% of the 
dynamical mass. Such large fracti ons of molecular ga s ma ss are 
also observed in local ULIRGs (ISanders & Mirabell 1 19961) as 
well as in distant SMGs (ITacconi et all 120061) . and indicate that 
the host galaxy of J0338+0021 is in an early evolutionary stage. 

4.2. The black hole-bulge mass ratio at z~5 

Based on the width of the Mgll 2798A line and the continuum 
intensity at A KSt = 3 000 A (Marinoni et al. in prep .), and by fol- 
lowing the prescription in lMcLure & Jarvisl (120021) . we estimate 
a black hole mass in this QSO of M B h ~ 2.5 x 1O 8 M0. We can 
infer an upper limit on the mass of a putative stellar bulge, by 
using the upper limit on the dynamical mass obtained above and 
subtracting the molecular gas mass (and assuming that the bulge 
is smaller than 3.5 kpc in radius). If the molecular gas disk is 



nearly edge-on, we derive MBH/Mb u i ge > 2.5 x 10 2 , which is 
substantially larger than the ratio observed locally, i.e. ~ 10~ 3 
dMarconi & Huntl l2003iPl . In order to have the lower limit on 
MfiH/Mbuige marginally consistent with the local value, the in- 
clination of the gas disk in J0338+00 21 must be about 20° , i.e. 
close to the average value found by ICarilli & Wangl (2006) for 
other QSOs with CO detection. However, as discussed above, 
the very broad CO emission of J0338+0021 relative to other 
QSOs suggests that the gaseous disk in the former is much more 
inclined. Moreover, one should keep in mind that the inferred 
MfiH/Mbuige is a lower limit. It is difficult to obtain more quan- 
titative constraints without higher resolution and higher sensi- 
tivity data. However, the current observations suggest that the 
MeH/Mbuige at high-z is higher than observed locally. This result 
is in agreement with the MBH/Mb u i ge ratio inferred for th e most 
distant QSO J 1 148+5251 at z=6.4 bv l Walter et all (120041) . Other 
indications of a hig her, with respect t o lo cal, MeH/Mbuiee mass 
ratio were found bv lPeng et al.l (120061) and lMcLureet al. (2006) 
in z>l AGNs. All of these results suggest that BH growth oc- 
curred on timescales shorter than bulge formation, and that the 
locally observed BH-bulge relation was achieved only at z<l. 

4.3. Star formation efficiency 

The far-I R emission is reg arded as a tracer of the star forma- 
tion rate dKennicuttl Il998l) . In QSOs the possible contamina- 
tion by dust heated by the AGN has been a hotly debated is- 
sue; however recent observations have shown that at least in the 
far-IR, the emission is gene rally due to star formation even in 
the case of powerful QS Os dSchweitzer et all l2006t iLutz et all 
120071: IWang et all 120071) . In J0338+0021 the observed L FIR = 
1.5 X 10 13 Lq implies a star formation rate of ~ 2500 yr" 1 , 
if usin g the LpiR-to-SFR conversion factor derived by lKennicuttl 
dl998h . The ratio Lfir/L^ is considered a measure of the star 
formation efficiency, since it is rela ted to the star formation 
rate p er unit of molecular gas mass (ISolomon & Vanden Bout. 
2005 ). Lfir/L^ is found to steadily increase with luminosity 
(e.g. Solomon et al.], 1 19971) . which is interpreted as an increas- 
ing star formation efficiency in the most powerful starburst sys- 
tems. J0338+0021 has a very high Lfir/L^ ratio, implying very 
high star formation efficiency. More specifically, in J0338+0021, 
the whole molecular gas content is expected to be converted 
into stars on a time scale of only ~ 10 7 yrs, i.e. a few times 
the dynamical timescale within the CO radius (td yn = R/v c < 
1 .5 x 10 7 yr). Cases like J0338+0021 are rare but not unique; in- 
deed similar "maximal starburst" systems are found among other 
hyper-luminous inf rared galaxies at low er redshift, both QSOs 
and starbursts (e.g. ITacconi et all l2006). 

4.4. Black hole accretion and molecular gas reservoir 

If we focus on type 1 (unobscured) AGNs, it is interesting to 
compare the molecular gas content, as traced by L™, with the 
optical luminosity, the latter being proportional to the Black 
Hole accretion rate. For type 1 AGNs with CO measurements 
(from Sanders et al l 1199 lb lEvans et all 120 06; Solomon et all 
1 19971: ISolomon & Vanden BoutL 12005b iMaiolino et all 11997). 
we have derived the rest-frame optical luminosity (/1L,i(4400A)) 
by using spectroscopic or photometric data available from the lit- 
erature at the observed wavelength closest to (l+z)4400A (gen- 

1 The value of MBH/Mbuige ~ 0.002 given in Ma rconi & HuntH2003h 
has to be lowered by a factor of 5/3 to account for more recent estimates 
of the bulge virial masses (Marconi priv. comm.). 
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Fig. 3. L' CQ versus (nuclear) AL A (440QA) for type 1 AGNs with 
CO measurements. Green triangles are sources at z<l, blue di- 
amonds are sources at l<z<5, red squares are sources at z>5. 
Large squares indicate the new measurements at z>5 reported in 
this paper. Crosses mark strongly lensed QSOs, for which dif- 
ferential magnification of the nuclear optical emission and CO 
emission in the host galaxy may occur. The black solid line in- 
dicates a linear fit to the data at z<l. Note the non-linearity of 
the relation and the likely flattening of the relationship at high 
redshift and high optical luminosities. 



erally near-IR data for QSOs at high redshift). For nearby Syl 
we considered only objects with nuclear measurements of the 
optical flux, to minimize the stellar light contamination. Fig. [3] 
shows L' co as a function of /1L / }(4400A), where sources in differ- 
ent redshift ranges are identified by different symbols and colors. 
Objects marked by a cross are strongly lensed QSOs: in these 
cases differential magnification may occur between nuclear op- 
tical emission and CO emission in the host galaxjQ. CO and op- 
tical luminosity clearly correlate, but the relation is non-linear. 
More specifically, by fitting the low-z data (z<l) alone, we ob- 
tain 



l; 



CO 



L Q K km s 1 pc 2 



= 51. 



AL A (440QA) 



0.72±0.08 



(1) 



High-z QSOs show indications for further flattening of the re- 
lation at high luminosities. In particular, QSOs at z>5 hint 
at a possible saturation of the CO luminosity at a few times 
10 10 L km s- 1 pc 2 . 

The correlation between BH accretion rate (/iL / |(4400A)) 
and molecular gas content in the host galaxy (L' CQ ) may be 
at the origin of the correlation between BH mass and stellar 
bulge mass observed in local galaxies (e.g. Ferrarese & Merrittl 



2 Note that high angular resolution observations rul e out significant 
lensing for all of the QSOs at z>5 consid ered here (Richard s et all 



lensing for all of the O S 
120061 12004 iFrev et all l2 



Prey et aul2005t IWhite et all 120051) . with exception of 
J 1630+40 12 for which no high angular resolution data are available. 



2000: [Marconi & Hunttl2003l) . The non-linearity of the relation 
is probably a consequence of the fact that, while BH accretion 
can be arbitrarily high (limited only by the Eddington luminos- 
ity), the molecular gas content is limited by the galaxy evolution- 
ary processes. In particular, the saturation of the CO luminosity 
in sources at z>5 may indicate that galaxy evolutionary mech- 
anisms cannot assemble more than a few times 10 10 M of gas 
in such early evolutionary phases of galaxy formation. Such re- 
sults are consistent with the finding of a saturati on also in terms 
of sta r formation rate in high-z, luminous QSOs (Mai olino et all 
120071) . 

Regardless of their interpretation, these results may explain 
the low galaxy-to-BH mass ratio observed in high-z QSOs (es- 
pecially at z=6.4), relative to local galaxies (Sect. l4~2i >. 



4.5. Low FIR/CO ratios are not common amongst high-z, 
luminous sources 

The detection of CO (5-4) and (1-0) in the far-IR weak 
(submm undetected ) radio galaxy TNJ0924-2201 at z=5.2 by 
iKlamer et alj (2005), suggested the possible existence of a sig- 
nificant population of high-z sources with large reservoirs of 
molecular gas but with little dust emission. These could be cases 
where dust had little time to form, or whose average dust temper- 
ature is extremely cold. However, the non-detection of CO in the 
two far-IR faint QSOs in our sample does not provide additional 
support for the existence of a large population of such objects, 
and TNJ0924-2201 remains the only case of exceptionally low 
Lfir/Lco (about a factor of 5 lower than sources with similar CO 
luminosity). As a consequence, strong far-IR emission seems to 
generally be a prerequisite for CO detection at high redshift. Of 
course the statistics are still extremely poor, and more observa- 
tions are required to investigate this. Moreover, we cannot rule 
out the possibility that in the two QSOs without CO detection 
the low-ionization UV lines provide a redshift which is offset 
by more than 1000 km s (which would move the CO line out 
of the old 3mm receiver band), although Mgll has a velocity 
generally consistent with the systemic ve l ocity within at most a 
few h undred km s _1 dVanden Berk et afl, 120011: iRichards et all 
2002), and our result on J0338+0021 supports this scenario. 

Acknowledgements. RM, PC, and MW acknowledge support from INAF. We 
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6 R. Maiolino et al.: Molecular gas in QSO host galaxies at z>5 

Table 1. Summary of physical properties and of the CO observations for the QSOs in our sample. 



Name 


ZMgii/cm] 


S v (850/im) 


S v ( 1.2mm) 


Lfir 


CO 


Ico 


zco 


FWHM 


I ' 


M(H 2 ) 






mjy 


mjy 


10 13 L o 


trans. 


Jy km s -1 




km s -1 


K km s pc 


10 10 M o 


J033829.31+002156.3 


5.027" 


11.9*2.0* 


3.7+0.3" 


1.5 b 


(5-4) 


0.73+0.09 


5.0267 


500 


2.7 x 10 1U 


2.2 x 10 1U 


J083643.5+005453.3 


5.774 c 


1.7±1.5* 


-0.4±1.0 S 


< 0.4 b 


(5-4) 


<0.41 






< 1.9 x 10 10 


< 1.5 x 10 10 


J163033.90+401209.6 


6.065 rf 


2.7±1.9" 


Q.8+0.6^ 


<0.8 e 


(6-5) 


<0.30 






< 1.0 x 10 10 


< 0.8 x 10 10 



Notes: upper limits are at 3cr; references for the redshift inferred from the low-ionization UV lines, for the 850/jm flux and and far-IR luminosity 
(corre cte d for our adop ted cosm ological constants) a re: " Marinoni et al. (in prep.). b iPriddev et all i2003h . c IStern et alj i2003h . rf lrwamuro etaD 
j2004l) . e lRobson et all j2004h . ^ iBertoldi et al] ( l2003l) . ^ iPetric et al.U2003l) , ICarilli et al.U200lh . 



